Dual Jets from Binary Black Holes 
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Supermassive black holes are found at the centers of most galaxies and their inspiral is a natural 
outcome when galaxies merge. The inspiral of these systems is of utmost astrophysical importance 
as prodigious producers of gravitational waves and in their possible role in energetic electromagnetic 
events. We study such binary black hole coalescence under the influence of an external magnetic 
field produced by the expected circumbinary disk surrounding them. Solving the Einstein equations 
to describe the spacetime and using the force-free approach for the electromagnetic fields and the 
tenuous plasma, we present numerical evidence for possible jets driven by these systems. Extending 
the process described by Blandford and Znajek for a single spinning black hole, the picture that 
emerges suggests the electromagnetic field extracts energy from the orbiting black holes, which 
ultimately merge and settle into the standard Blandford-Znajek scenario. Emissions along dual and 
single jets would be expected that could be observable to large distances. 



Introduction: Among the most spectacular of astro- 
physical events are observations of tremendous amounts 
of energy careening out along highly collimated jets from 
a localized central region. Within this central region is 
the engine behind the jets, generally believed to be a 
spinning black hole which helps to convert binding and 
rotational energy into kinetic and thermal energy of the 
surrounding plasma (see e.g. [1 ). This basic picture is 
widely accepted, but a detailed understanding of these 
systems remains elusive in part due to the inability of 
detecting clean electromagnetic signals from the depths 
of the central engine. Such understanding is particularly 
needed to understand phenomena such as gamma ray 
bursts and AGNs, which radiate a tremendous amount 
of energy and, in some cases, are thought to play a fun- 
damental role in the formation and evolution of massive 
galaxies |2]. 

In addition to intense efforts in various electromagnetic 
bands, observations in the gravitational wave band will 
soon be added. Because gravitational waves propagate 
essentially unscattered through the universe and are most 
strongly produced in highly dynamical regions of strong 
spacetime curvature, their observation promises to open 
up a new complementary window to the universe, allow- 
ing for multi- messenger astronomy ^ |4^. Astrophysi- 
cal systems expected to be bright in both gravitational 
and electromagnetic spectrums are therefore attracting 
intense attention. 

A particularly interesting scenario involves the merger 
of a supermassive binary black system which could be 
observable with LISA to redshifts up to z = 5 and be- 
yond [5 . Such a system is expected as a natural conse- 
quence of galaxy mergers, whose individual black holes 
eventually settle into an orbit sufficiently tight that the 
binary's subsequent dynamics is governed by the grav- 



itational radiation timescale and disconnected from the 
properties of the disk [6]. As they orbit, they emit grav- 
itational waves which carry off both energy and angular 
momentum, driving the black holes to merger producing 
bright gravitational waves. The expected gravitational 
wave signals from such a system, their impact on gravita- 
tional searches and their exploitation for physics mining 
through data analysis are now well understood thanks 
to concentrated efforts on the analytical and numerical 
fronts (see e.g. [7 and references cited therein). Such ef- 
forts have studied the dynamics of orbiting binary black 
holes in vacuum, which is a good approximation as the 
black hole's inertia is orders of magnitude larger than any 
other one. 

A logical next step is to comprehend possible electro- 
magnetic counterparts which would allow for studying 
key systems via both bands. Such a possibility is quite 
strong in the context of supermassive binary black hole 
mergers, because in the merger process the black holes 
will interact with at least: (i) a circumbinary accretion 
disk [6 , (ii) remnant gas in between the black holes ^[9], 
(iii) a magnetosphere [10] and give rise to possible elec- 
tromagnetic emissions (e.g. [TTHT7]: pHSOj: [2Tti23] ^. 
As discussed in [6 , as the binary tightens, a common 
circumbinary disk will be formed within which the black 
holes will eventually merge. Such a disk will typically be 
magnetized, anchoring a magnetic field which will perme- 
ate the region inside the disk containing the black holes. 
Once the black holes merge, the generic outcome will 
be a spinning black hole within the magnetic field an- 
chored by the disk. It is precisely this system studied 
in the seminal work of Blandford and Znajek [10 and a 
large body of subsequent works. In this model, the spin- 
ning black hole immersed in an external magnetic field 
accelerates stray electrons to energies sufficient to pro- 



duce cascading pair-production which suppUes the envi- 
ronment with a plasma [1]. Due to the neghgible fluid 
inertia, the bulk of such plasma moves freely without ex- 
periencing any Lorentz force, and so the system can be 
cleanly studied with the so-called force-free approxima- 
tion [TOl [24] . This approximation, together with the as- 
sumption of stationarity and axisymmetry, allows one to 
determine that a net outward electromagnetic flux is pro- 
duced, which extracts rotational energy from the black 
hole. Furthermore, the electromagnetic energy flux is 
highly collimated and eventually will be transferred into 
kinetic energy of the plasma, accelerating charges and 
producing synchrotron radiation. Support for this basic 
picture has been provided in recent years through numer- 
ical models [25ti29]. 

The Blandford-Znajek model outlined above is not 
directly applicable to the very dynamical merger of two 
black holes, thus the interaction of electromagnetic fields 
and plasma and their possible emissions in a such system 
has remained a mystery. Our studies indicate that the 
dynamical behavior of the system: induces a collima- 
tion around each intervening black hole -generating a 
toroidal magnetic field — ; amplifies the electromagnetic 
field strength (when comparing to the single black 
hole case), and emits a strong electromagnetic burst 
associated with the merger. The collimated area can 
thus become a channel for accelerated particles which, 
in turn, can emit observable electromagnetic radiation, 
in particular of synchrotron type. Furthermore the 
behavior of the EM fields, tightly tied to the binary's 
dynamics, make them strong tracers of the spacetime 
and the collimation tubes can be seen as a cover for the 
familiar "pair of pants" description of the black holes' 
event horizons [30] , 

Implementation details: We study the black hole 
binary via numerical simulations which implement the 
(coupled) Einstein-Maxwell system of equations. Also 
coupled to the system is a low density plasma for which 
the inertia is negligible compared with the electromag- 
netic stresses. The plasma essentially experiences no 
Lorentz force and so the electromagnetic fields are well 
described by the force free condition Pem^+Jem A5 = 
(with pem, Jem the charge and current densities). The 
role of the plasma is then to provide charge and cur- 
rent densities. The complete system is implemented 
using finite difference approximations within a compu- 
tational infrastructure that provides distributed adap- 
tive mesh refinement [31 with seven levels of refine- 
ment. Our computational domain is defined by x* G 
[-1207 Ms AU, 1207 Ms AU] with a finest resolution of 
A = OAMsAU, where M is the total mass of the system 
in terms of the solar mass Mq and Ms = M/10^Mq. 
We adopt fourth order accurate spatial discretizations 
and third order Runge-Kutta integration in time, while 
the singularity inside each black hole is excised from the 



computational domain. Since the region inside the black 
hole is causally disconnected from the outside, this pro- 
cedure does not affect the physics obtained there. More 
details of our techniques for implementing Einstein equa- 
tions can be found in e.g. [2T|[32H34] and for our adoption 
of the force free condition we follow closely the approach 
presented in [25 . 

The circumbinary disk is assumed to lie outside the 
computational domain and so its details are unimpor- 
tant. Such an assumption is strongly justified by the fact 
that the viscosity in the disk is such that the disk can- 
not keep pace with the shrinking orbit of the black holes, 
and as a result the disk "freezes-in" at a radius typically 
considerably larger than our computational domain [6^. 

The disk's role is therefore to anchor the magnetic 
field which is incorporated as boundary and initial 
conditions describing a dipolar magnetic field with 
strength Bq in the region around the black holes. We 
adopt black holes in a quasi-circular orbit (which is a 
reasonable approximation as gravitational wave emission 
circularizes the orbit). We study late orbiting stages and 
merger, and so we concentrate on an initial separation 
that gives rise to over one orbit before the merger takes 
place. This scenario allows us to study the transition of 
the binary from orbit to merger into a final black hole. 
Furthermore, we adopt equal-mass {Mbh = M/2, with 
a radius Rbh = lOMsAU) non-spinning black holes 
to disentangle orbitally-driven effects from those that 
would be induced by any individual spins (i.e. individual 
BZ effects that would be present around each spinning 
black hole). 

Results: As mentioned, gravitational waves from this 
system, in which the plasma and electromagnetic inertia 
are orders of magnitude smaller than that of the black 
holes, are well understood by now (see [7 and references 
cited there in) and our results reproduce the expected be- 
havior. The electromagnetic field behavior on the other 
hand, is completely new and we concentrate on it here. 
The global dynamical behavior is illustrated in Fig 1. 
The two black holes orbit and merge within about an or- 
bital period {^ 13.6 Mg hours), after which the black fi- 
nal hole region radiates its excess structure, settling into 
a Kerr black hole with spin a = J c/{G M'^) c^ 0.67. 
To study the system in more detail, we monitor differ- 
ent quantities: The Newman-Penrose scalars ^2 and ^4 
(the square of the former and the integral of the lat- 
ter describe the electromagnetic and gravitational en- 
ergy fluxes respectively); the Poynting flux; the function 
Qf = Fro/Frcf) (which in axisymmetric, stationary sce- 
narios describes the angular velocity of magnetic field 
lines [10]) and the electric and magnetic field topologies. 
To describe the behavior of the system, it is easier to con- 
centrate on the different stages identified in the binary 
black hole dynamics: 

In the first stage, before the merger takes place, the 










(a) -11.0 Ms hrs 



(b) -3.0 Ms hrs 




(a) -8.2 Ms hrs 



(b) 4.6 Ms hrs 



FIG. 2: Some representative electric and magnetic field lines, 
together with the electromagnetic rotation frequency Qf- 



(c) 4.6 Ms hrs 



(d) 6.8 Ms hrs 



FIG. 1: Electromagnetic energy flux at different times. The 
collimated part is formed by two tubes orbiting around each 
other following the motion of the black holes. A strong 
isotropic emission occurs at the time of merger, followed by a 
single collimated tube as described by the Blandford-Znajek 



black holes act as stirrers for the surrounding plasma. 
Their orbital dynamics and strong curvature affect the 
electromagnetic field in the close neighborhood of each 
black hole, inducing both a poloidal electric field and 
a toroidal magnetic field -both scaling as vBq- (as in 
models for magnetospheric interactions of binary pul- 
sars [35l [36]). Fig. 2 illustrates field lines correspond- 
ing to E and B as well as the VLp. The induced time- 
dependent topology gives rise to a net Poynting flux 
aligned/antialigned with the orbital momentum vector 
around each black hole, with an EM frequency given by 
^F ^ ^orb/5. As a result, despite the fact that the black 
holes are not spinning, there is a strongly collimated elec- 
tromagnetic flux of energy dominated by an tti = 2 mul- 
tipolar structure with a time-dependence determined by 
the orbital motion. 

During the merger, which lasts about ~ TMg hours 
-from late orbiting, through the plunge phase to the for- 
mation of a highly distorted black hole- the highly non- 
linear dynamics translates into a significant increase in 
both electromagnetic and gravitational energy radiated. 
As expected, a common collimated tube arises and a nat- 
ural transition from tti = 2 ^ is observed. In addition 
to the collimated flux of energy, a rather isotropic flux is 
also emitted whose energy is much smaller than the col- 
limated one during the inspiral but of the same order at 
the merger, where both these energy fluxes reach a peak 



and the collimated part doubles in magnitude. 

Afterwards, the late stage is described by a single black 
hole which, after a relatively short time, is well approxi- 
mated by a Kerr black hole. The remnant black hole ra- 
diates gravitational radiation described by quasinormal 
modes which decay to zero exponentially as the black 
hole settles into a Kerr configuration. Thus, the post- 
merger behavior of the electromagnetic field behavior is 
increasingly better represented by the Blandford-Znajek 
process for a spinning black hole with a :^ 0.67. As a re- 
sult, the collimated electromagnetic energy flux does not 
decay to zero, rather it approaches the value predicted by 
the Blandford-Znajek model. For a single spinning black 
hole, this energy flux evaluated at the horizon goes like 
Fem ^ Rbh^f(S^h — ^f)B^, where Qh is the rotation 
frequency of the black hole (which is similar to the orbital 
velocity at the merger). In this case, we numerically find 
that the EM rotation frequency for a single black hole 
relaxes to Qf ^ ^h/2. 

Energetically, the system radiates gravitational waves 
primarily through I = \m\ = 2 modes. These waves 
display a chirping behavior as the orbit tightens, fol- 
lowed by exponential decay after the merger (see e.g. [7]). 
Overall, the system radiates c^ 2.5% and ^ 16% of the 
rest mass energy and angular momentum (at the initial 
separation) respectively. In the electromagnetic band 
the radiation profile displays a more complex structure. 
Throughout the orbiting stages, the electromagnetic ra- 
diation exhibits a strongly collimated character along the 
magnetic field lines in the region around the individual 
black holes, together with a weaker isotropic emission. 
These features are illustrated in Fig. 1 in which the flux 
of electromagnetic energy is shown at four different times 
during the evolution. In the early stages the black holes 
stir the surrounding plasma, leading to a clear collima- 
tion of electromagnetic fiux induced by the orbiting black 
holes pulling the EM field lines. These collimated tubes 
twist about each other as the black holes proceed. When 
the merger takes place, the collimated tubes merge into 
one and acquire a rather smooth structure around the fi- 
nal black hole's ergosphere. In addition to the collimated 
energy fiux, a strong burst of isotropic electromagnetic 
radiation is produced at the merger. 






(a) -8.2 Ms hrs (b) -5.5 Ms hrs (c) -3.0 Ms hrs 




(d) 2.0 Ms hrs 
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(e) 4.6 Ms hrs 




(f) 6.8 Ms hrs 



FIG. 3: Projection of the radial electromagnetic flux on a 

sphere located at r = 20M. The structure of the tubes and 

the transition from a?7i = 2toa?7i = mode are clearly 

displayed. 
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FIG. 4: (Left) The gravitational and the (collimated) elec- 
tromagnetic luminosities for the case Bo = 10^ G and M = 
10^ M0. The gravitational luminosity increases more than an 
order of magnitude at the merger and then decays to zero. 
The electromagnetic luminosity during the inspiral is only a 
factor ^ 2 smaller than at the merger. After the peak, it 
decays to a constant value given by the Blandford-Znajek 
mechanism for a single spinning black hole. (Right) Elec- 
tromagnetic (yellow-red) and gravitational (green-blue) radi- 
ation from the system at 4.6Ms hrs. 



Figure 3 illustrates the behavior of the electromagnetic 
energy flux over a sphere located at r = 196 MsAU. 
A clear collimation is observed which is evident in the 
"bright-spots". These spots are symmetric with respect 
to the orbital plane and revolve around each other as a 
result of the orbiting until that they merge into a single 
one along the poles. In addition to these collimated flows, 
the figure also illustrates the energy flux through differ- 
ent latitudes and the transition from the m = 2 structure 
to m = 0. 

Figure 4 provides the total energy flux in both the elec- 
tromagnetic and gravitational wavebands for a particu- 
lar astrophysically relevant case with M = 10^ Mq and 
B = 10^ G (i.e. not exceeding the Eddington magnetic 
field strength B c^ 6 ■ 10^ M~^^^ G [37 ). 

A clear sweep upwards is apparent in both channels 



as the merger takes place. Afterwards, both decrease 
rapidly. However it is interesting to note that while the 
gravitational flux essentially shuts off in a short time (as 
the system relaxes to a stationary black hole), the elec- 
tromagnetic flux rises up to the level predicted by the 
Blandford-Znajek mechanism. 

Overall, the total (Poynting flux) radiated elec- 
tromagnetic energy behaves like E^^^/{Mc^) = 
10-^"^ [Msf [B^f (with ^4 = B/{l{)^G). This is about 
an order of magnitude larger than the same system stud- 
ied in the electrovacuum case (i.e. not including the 
plasma) [211 122] ^ indicating that the plasma taps the or- 
bital/rotational energy from the system more efficiently. 
As studied in [38 , the isotropic radiation transport of 
magnetic field perturbations can be transmitted to the 
disk in the form of magnetosonic waves when the ratio 
between the viscous time scales of the disk and the trans- 
port magnetic timescale is above :^ 1. Thus, in [38] it was 
indicated that an electromagnetic counterpart to a grav- 
itational wave event could be determined by a variability 
induced in a possible accretion mechanism taking place. 

The results highlighted in the current work indicate 
that a more important observational prospect is achieved 
by the jets induced by the binary's dynamics. These 
jets have a (Poynting flux) luminosity during the last 
hours before the merger of 2 — 4 • 10^^ erg/ s ~ 0.0015 — 
0.0030 L^;^^^ (with LEdd the Eddington luminosity). This 
electromagnetic energy can be transferred to kinetic en- 
ergy of the plasma, which will radiate through syn- 
chrotron processes in frequencies around the GHz region. 
Such a flux of energy implies it would be possible to ob- 
serve these systems to z c^l^ with missions such as IXO 
and EXIST ^ depending on the efficiency of the photon 
emission process. As Fig 4. illustrates, the system will 
be sufficiently bright for hours before the merger and 
remain so through merger. This radiated power has a 
time dependence given by (i^orb ^orb)^^^ and the flux 
exhibits a clear transition tti = 2 ^ at the time of the 
merger. Thus, prospects for detecting pre-merger and 
prompt electromagnetic counterparts are certainly inter- 
esting. Furthermore, joint detections in both electric and 
gravitational wave bands (the power of the latter scales 
as R^^^Vt^^^^M'^) are therefore quite probable as LISA 
will be capable of observing supermassive binary black 
hole systems for weeks to months before merger and con- 
siderable time afterwards. 

Conclusions: Observational prospects of detecting 
gravitational waves with LISA from supermassive black 
hole mergers will be excellent [39 . Indeed, EISA's superb 
noise spectra will allow for clean determination of gravi- 
tational waves and extraction of key physical parameters 
of the system such as certain combinations of the their 
masses and spins, orientation, position in the sky and lu- 
minosity distance. Gravitational waves from supermas- 
sive binary systems could be detected as far as redshifts of 
2: = 5 — 10 [5]. However, EISA's localization in the sky-an 



ellipse with a few arcminutes to a few degrees in size- will 
not be sufficient to unequivocally determine the source's 
position [39l |40] , and it is here where an electromagnetic 
counterpart would aid tremendously in determining the 
location. Precise localization and confrontation of the 
observed signals will allow for tremendous physics pay 
off: general relativity, cosmology, astrophysics and even 
efforts towards a quantum theory of gravity will likely 
enjoy a revolution in their understanding J39l |4T| [42] . 

As we have shown here, binary black hole mergers 
could give rise to strong electromagnetic output with 
(Poynting flux) luminosities as high as a few 10^^ ergs/s. 
This would correspond to an isotropic bolometric flux of 
Fx ^ 10"-^^ erg/ (cm^ s). Such a flux could be detected 
to a redshift of z = 1 and possibly larger if accounting for 
anisotropics on the one hand, but on the other hand de- 
pending on the photon emission process efficiency on the 
other. Interestingly, our work indicates that a marked 
transition in the electromagnetic flux from m = 2 ^ 
will take place around the merger time which will leave 
its imprint in possible observable signals. Furthermore, 
we have shown that even non-spinning black holes could 
give rise to interesting levels of output, which becomes 
even stronger as the merger takes place. Certainly, higher 
power would be expected if any of the individual black 
holes is spinning and/or the magnetic field strength is 
larger in the region of the black holes either as sourced 
by the disk or further enhanced by the merger process. 
In the latter case, a tantalizing option would be for such 
enhancement to be driven by an analog of the BZ effect; 
i.e. by significantly tapping orbital rotational energy be- 
fore the merger takes place as in the single black hole case 
where rotational energy from the black hole is extracted. 
A requirement for this scenario would be that the black 
holes orbit within an ergo-region and before the plunge 
so that there is sufficient time to efficiently tap into this 
source of energy. As argued in [22 , this scenario appears 
unlikely unless the black holes are significantly spinning 
and aligned with the orbital angular momentum; oth- 
erwise, the black holes would be in the plunging phase 
when crossing (an estimate of) the ergosphere. While 
the magnitude of spins in supermassive black holes is 
still largely unknown -gravitational wave detections will 
be the ultimate tool to understand this issue- alignment 
of spins with the orbital angular momentum is expected 
in a large class of systems [43 . 
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